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INTRODUCTION 



The isotropic multichannel spin-S* Kondo model de- 
scribes the XXX-\ikc spin scattering between non- 
interacting spin-1/2 fermions in m channels and a sin- 
gle localized spin-S* impurity. It has been proposed by 
Nozieres and Blandir*i to account for the orbital struc- 
ture of spin-S* impurities in metals. Physical realizations 
are discussed in the review articla^i. An exact solution 
to this model was found by Tsvelick and WiegmannMii 
using the Bethe Ansatz (BA). By applying the thermo- 
dynamic Bethc Ansatz (TBA), the thermodynamics has 
been calculated by Tsvelicl*& and Andrei and Destri^. 
In the TBA approach, the free energy is encoded by 
infinitely many coupled non-linear integral equations. 
Recently, SchlottmannS*^ and Zarand et aliS obtained 
the free energy for the anisotropic multichannel spin-S* 
Kondo model with XXZ-\ike exchange in the TBA ap- 
proach. Schlottmann observed a quantum critical point 
for the overscreened model in the limit of low temper- 
atures, Zarand et al focused on the underscreened two- 
channel model. Here, we present a novel exact solution 
in which the free energy is encoded in a finite set of 
non-linear integral equations (NLIE), including both the 
anisotropic and isotropic models for arbitrary S and to. 
On the one hand, we confirm known results by this novel 
solution, on the other hand, we obtain farther reaching 
results as far as ratios of high- and low-temperature scales 
are concerned. 

In a recent publicationii, we proposed a lattice path 
integral approach to an Anderson-like impurity in a cor- 
related host. We obtained this model from the Hamil- 
tonian limit of a gl(2|l) symmetric transfer matrix. The 
free energy contributions of both the host and the im- 
purity were calculated exactly by combining the Trotter 
Suzuki mapping with the quantum transfer matrix tech- 
nique. We showed that the model allows for the Kondo 
limit, i.e. a localized magnetic impurity in a free host 
with linearized energy-momentum relation. In the Kondo 
limit, the symmetry of the coupling between host and im- 
purity is reduced from gl(2|l) to su(2). In the following, 
we exploit this observation to propose a UqSu(2) sym- 



metric quantum transfer matrix (QTM), whose largest 
eigenvalue yields the free energy contribution of the im- 
purity in the thermodynamic limit, in close analogy toii. 
The g-deformation accounts for an XXZ-\ikc spin ex- 
change between the localized moment and the host, in- 
cluding the isotropic limit q ^ 1. The transfer matrix 
is built from local i?-opcrators, each acting in the ten- 
sor product space of two UqSu(2) modules, carrying an 
(?7i+l)-dimensional irrep and an (Z + l)-dimensional irrep 
of UqSu(2), respectively. Here m denotes the number of 
the electronic host channels and Z = 25" is twice the im- 
purity spin. By making use of analyticity arguments, the 
free energy is represented by (max[Z, m] + l)-many NLIE. 
These equations arc investigated both analytically and 
numerically for the accurate calculation of scales in the 
limits of low and high temperatures. 

This article is organized as follows. In the next section, 
the free energy contribution of the impurity is determined 
as the largest eigenvalue of a certain QTM. The third and 
fourth sections are devoted to the calculation of thermo- 
dynamic equilibrium functions and scales in the limits 
of low and high temperatures, respectively. In the last 
section, we summarize the results. 

In all that follows we set fcs = 1, and g/Js ~ 1, where 
kB is Boltzmann's constant, g is the gyromagnetic factor 
and fj,B is the Bohr magneton. An index i (h) denotes 
quantities pertaining to the impurity (host). 



2. CALCULATION OF THE FREE ENERGY 

Let Vq ^'' be the module carrying the l-dimensional ir- 
rep of UqSu(2); in the limit q = 1, V^''^ carries the 
Z-dimensional irrep of su(2). Consider the matrices 

Rf'^\x) eEnd^vJ'' ® fJ"^) constructed such that the 
Yang-Baxter-Equation (YBE) is fulfilled: 
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The R'i'"'\x) are obtained by fusing a lattice of m x I the trace over the auxihary space of Tg'' 
many operators Rq^'^'^ . Then Vq'^^ is the subspace of com- 



In view of 



pletely symmetric tensors in 



Va 



(1) 



with k = /, m. 



Ini^, the exphcit expression of Rq'™\u) is given. If 
in the isotropic single channel case, the spectral pa- 
rameter is chosen to take the special value xq := 



i+i 



- tan J {I + 1) , the matrix R^^'^^ (xq) simplifies toi^ 



R['''\^o) 



exp [12 JS ■ 



(1) 



The logarithm of equation is the spin exchange oper- 
ator between the impurity spin S = 1/2 and one electron, 
where the spin-exchange coupling J is parameterized by 

Xq. 

It has been argued by Tsvelick et al^ that the exchange 
in the multichannel Kondo model is equivalent to an ex- 
change model for spin-7Ti/2 electrons scattered by the 
spin-S' impurity. This equivalence has been proven by 
the equality of the exact solutions of both models at zero 
temperature. We expect that more insight into this ques- 
tion is possible by a lattice path integral formulation in 
analogy to the 5* = 1/2, m = 1, 7 = caseii. We leave 
this task as a future challenge. Here, we rely on the above 
mentioned equivalence, so that the Hamiltonian density 
is given by: 



H = Hh + Hsd + 'Hex 
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the definition of the Hamiltonian in equation Ij2b|) . the 
auxiliary space a is identified with the impurity space and 
the quantum spaces 1, . . . , iV correspond to the electrons. 
We will now show that under the substitution N — > \(3D, 

the largest eigenvalue K^i'^\xq) of t[^'^\xq) is directly 
related to the free energy contribution of the impurity fi 
found inii, 
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Here D has the meaning of a bandwidth parameter of the 
host. A magnetic field is introduced by twisted boundary 
conditions. 

The largest eigenvalue in equation is expressed 
through auxiliary functions b, b (*B = 1 -I- b, *8 = 1 -I- b) 
in analogy toi4. After the substitution N if3D, these 
functions are given by the non-linear integral equation 
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e d/c and k,±[x) := K(a;± 



-00 2 cosh ■ 

m). An equation similar to ^ holds with b ^ b, ^ —h 
and K_ K+. The free energy is obtained from the 
auxiliary functions by 
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These equations are equivalent to the corresponding TEA 
equationsi^, as can be shown by fusion techniques de- 
scribed ini^. Proceeding as ini^, the universal limit 
D — !■ (X) is taken by defining the temperature scale 



Tk := 2De~ 



This results in the Z3-indcpendent equations 



(6) 



where xq and q parameterize the coupling constants. In 
the following, q is given by a real constant 7, q =: exp(i7), 
such that 7^0 induces an Ising-like anisotropy, Jj^ < 
J||. The matrix Rq'"^^ acts in the electronic space with 
operator-valued entries in the impurity space. Only the 
totally symmetric combination of spinors, with total spin 
m/2, interacts with the impurity, so that the spin indices 
in equation (|2bp are u, a' = —m/2, . . . , m/2. 
Consider the monodromy matrix 



dh 

Inb(x) = -e^ + ^-h [K*ln<B-K_*lnS] (x) 
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in complete agreement with results fromii. Thus we ar- 
gue that the impurity contribution to the free energy 
of the anisotropic multichannel spin-S" Kondo model is 
given by 
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In the last line the transfer matrix r, 
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acts in the tensor product Va'' ®Vn"^\ 



is defined as 
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where the universal limit Z) ^ 00 is implicit. 

The eigenvalue Ag'™'' is found by standard Bethe- 
Ansatz technique9iSii&. FoUowingil, one introduces suit- 
ably chosen auxiliary functions in order to exploit the 
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analyticity properties of the largest eigenvalue. We re- 
strict ourselves to the anisotropy parameter 



T <^ 2^ ■ Tmax ; 



(8) 



where n := max(Z,m). We define a temperature scale 



from which the isotropic case © is recovered by scaling 
xq 7x0. One is left with the following system of NLIE 
for the auxiliary functions yj = 1+Yj, (j = 1, . . . , ?i — 1), 
b„ = 1 + *B„ and b„ = 1 + 
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and the free energy contribution of the impurity is given 

by 



-T ■ 



roc Inyi(x) 

27rcosh(2-+ln ^) 



InYijx) 
tsh(^a;+ln 

27r cosh( x-{-\n tS— ) 



dx , I < m 



■ dx , I > 



.(14) 



Again, the magnetic field is introduced by twisted bound- 
ary conditions. The integration kernel k depends on n 
if 7 7^ 0. The asymmetric driving term — e^ in the m-th 
equation gives rise to different asymptotes in the limits 
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and similarly for b„ with h ^ —h. 

The above presented derivation of the closed set of 
finitely many integral equations for the impurity's free 
energy is speculative in the sense that an observation 
made for the well-understood case S = 1/2, to = 1 within 
a path integral approachii is extended to the case of gen- 
eral S and m. We intend to report on the generalization 
o£^ in a future publication. Our claim, however, that 
(|10|l - p3|) . (|14|) hold also in the general case follows by 
a simpler reasoning. By use of the algebrai^*ii of fused 
transfer matrices one can show that our key equations 
are equivalent to the infinitely many TBA equations o& 
in the isotropic case and to the equations o&S in the 
anisotropic case. As an illustration, figures [ll^ show the 
entropy and specific heat for some special cases. The 
curves for 7 = agree nicely with similar graphs in2ii^; 
those for 7 ^ are novel. 

We expect that in the framework of a generalized lat- 
tice path integral approach, additional electron-electron 
interactions in the host must be introduced to keep the 
model integrable. These interactions may lead to both 
a renormalization of the Fermi velocity Vf and the elec- 
tronic g-factor, analogous to the known exact solutionis. 
These additional interactions will be the concern of fu- 
ture work; the analysis in this work is based on the model 



4 





- 








— m=\ 




m=2 




-— m=3 




m=4 




- - m=5 




— Asymptotes 



10 




FIG. 1: Entropy for impurity spin S = 5/2 and differ- 
ent ciiannel numbers. At low temperatures, tfie asymptotes 
ln(m + 1) depend on the diannel number, at high tempera- 
tures, the common asymptote ln(2S -I- f ) is approached. 



3. LOW TEMPERATURE EVALUATION 

3.1. 2S = m: Exact screening 

The first non- vanishing order in a low temperature, low 
field expansion of the free energy in the exactly screened 
case can be obtained by expanding the corresponding 
integral in the following form 
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The temperature and the field are supposed to be small 
compared to Tx- As can be seen from equations (|15|l . the 
integral in equation H16|l only exists for I ~ m. The above 
approximation makes sense if T, /i Tk- The integral in 
equation H16|l can be done exactly, by a generalization of 
a method used iniS leading to dilogarithms and applied 
for instance to the spin-S" Heisenberg chain irkiiSS. 
In the notation of equation H1U|I . consider the integral 
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From the symmetry of the integration kernels one de- 
duces 
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On the other hand, the integral /(T, h) is calculated as 
in^ by using dilogarithmic identities, the result is 



m + 2 ^ ' 

Combining equations p8|l . 1)19(1 . one gets the first term 
in an expansion of the free energy for low fields and tem- 
peratures: 

MT,h) : 



2Tk m + 2 4(7r - mj) Tk 



7^- (20) 



FIG. 2: Specific heat for impurity spin S = 5/2 and different 
ciiannel numbers. The inset shows the Fermi liquid behavior 
limr^o C{T)/T = const, for 25 = m = 5. 



From equation H2U|I . 
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These relations exhibit clearly the role of Tk as a "low 
temperature scale" . According to equation 1)21(1 , the im- 
purity contribution to the specific heat and magnetic sus- 
ceptibility are Fermi liquid like at T, h <^ Tk in the ex- 
actly screened case I ~ m. This is the regime of "strong 
coupling": The anti-ferromagnetic spin exchange leads 
to the formation of a many particle state between the 
impurity and the host electrons, which screens the mag- 
netic moment of the impurity. Elementary excitations 
of this bound state are Fermi like. Nozieresi^ built up 
a phenomenological Fermi liquid theory to describe this 
regime. This Fermi liquid behavior is to be compared 
with the host. It consists of spin-l/2-fermions of m non- 
interacting channels (or of m fiavors), so that the density 
of states is enhanced by a factor of m. 
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The coefficient of the linear T-dependence of Ci [Gh) is 
denoted by 5i {5h)- The low-temperature Wilson ratio R 
is defined and calculated as 
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The lower bound 2 is reached for m ~ 1, 7 = 0. The 
striking feature in comparing equations 1(21(1 . 1(22(1 is that 
Ci is reduced by a factor of 3/ {2m + 4) in comparison to 
C/i, if the constant ph is chosen such that Xi = Xh for 
7 = 0. This may be interpreted by the localization of 
the impurity: Contrary to the host electrons, it does not 
move, so that the specific heat is reduced. 



3.2. Linearization 

In this section, the ground-state and the lowest T- 
dependent contribution to the free energy are calculated 
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FIG. 3: Entropy for m = 5 channels and 

spins 5* — 1/2, ...5/2. In the low tempera- 

ture limit, the entropy takes the unusual value 
ln{sin [71(25"+ l)/(m + 2)] /sin7r/(m + 2)}, equation (ITFEl . 
whereas for high temperatures, the asymptotes are given by 
ln(2S+ 1). 



by linearizing the NLIE for h ^ 0. From equations H15|) 
one observes that for h =/= and f3 —^ oo in the limit 
X —oo the auxiliary functions scale with Ph. Es- 
pecially, lim/32>i In 05; = 0-1-0 (e~'"'), so that one can 
neglect In 03; with exponential accuracy. 
One introduces the scaling functions e^. 



lnyj{x) =: PhKm ej{-x + \n(3hK^) j = 1, 
\nbi{x) =: l3hKmei{-x + In phK„r) , 



where Km is defined by lny,„(ln/3/iiir,„) = 0. The shift 
in the spectral parameter is performed in order to deal 
with functions which have a zero in the origin. 
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so that one linearizes the logarithms 
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(ihKm^j + C'(exp(— const. /3/i)), j > m 

The linearization is done with exponential accuracy for 
£j>rm but with only algebraic accuracy for ej<m because 
of the constant asymptotes of the latter functions. The 
first T-contribution for these is calculated below within 
a different linearization scheme. Here, the lowest T- 
contribution for the exactly and under-screened single 
channel case ^ > m = 1 is determined. For 7 = 0, the 
results o£i^ are confirmed. 

It is convenient to define a matrix A^^ 
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FIG. 4: Specific heat for m — 5 channels and spins S 
1/2,... 5/2. 



By inserting the e-functions into the original set of NLIE 
and using equations (|26|) . H27|l one obtains 
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where we have defined the driving term 
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The terms ej^o contain the ground-state, and from the 
rest, the lowest T-dependent contributions are obtained. 
We will calculate the latter explicitly for / > m = 1. We 
made use of the shorthand-notations 

e,{x) e{x) =: e+{x) , e,{x) e{^x) =: ej {x) 

for functions in direct space. Their Fourier transforms 
are denoted as 
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hence an index + denotes analyticity in the upper half of 
the complex fc-plane. Define a new energy scale 
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Then the free energy, magnetization and specific heat are 
given by 
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FIG. 5: The entropy for the under-screened case S = 
1, m — 1 with fields h — 0,0.1, 1, 10 and anisotropies 
7 = 0, 0.97max — 0.707 (for 7max, see equation JHJ). 




FIG. 6: The specific heat for the under-screened case S = 
1, m = 1 with fields h = 0, 0.1, 1, 10 for anisotropies 7 = 
0, 0.97max — 0.707 (for 7max, see equation JHJ). 



Obviously, the Wilson ratio Ci/{Txi) is a function of 
h alone. In the following, the magnetization and the 
susceptibility are calculated from e;.o, and the specific 
heat from Ui, defined in equation 129|l . 



3.3. Under-screened and exactly screened cases, 

2S> m 



The linear system H28|) is solved by inverting the matrix 
A~^{k). This gives the matrix A{k) with elements 
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where Uij is the adjunct to 

quantities are seen to satisfy recursion relations, which 
can be solved explicitly with the results 
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Other matrix elements are not needed. Then 
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Eq. H36|) is solved by factorizing Am^m into two functions 
G±, G+ (G_) being analytic in the upper (lower) half of 



the complex /c-plane (Wicner-Hopf factorization). 



The result is 
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where k is the Fourier transform of the kernel defined in 
equation 112|) for / = 1 and 
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The constant a is chosen to be a 
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Y ln(l — tt/ (7771)). From equation (|36|l the constant 
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With this information and equations (|30|l we know 
the magnetic field scale Th whose role is revealed in var- 
ious expansions of the magnetization m{h) below. 



Ground state 

The ground state contribution to the free energy equa- 
tion (EH) is given by inserting ei^ from equation H37|) with 
equation H39|) . 

The magnetization is written down from equation H32|l . 
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The integral can be caleulated by closing the contour 
in the lower (h > Th) or in the upper {h < Th) half 
plane. This results in two power series in T^/h for h > Th 
and in h/Th for h < Th with integer and non-integer 
powers, depending on the poles of the integrand. These 
calculations are straightforward but lengthy. Especially, 
in the limits h Th, h ^ Th, one extracts the constant 
values 
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Thus for h ~ 0, a non-integer rest-spin remains in the 
case 7 ^ 0, Z > m and ^.-dependent corrections are of 
non-integer powers of h. The constant terms in H43|l 
agree precisely with the asymptotes, equations (I15II and 
witl&S, where the same model was investigated by TBA- 
techniques in the limits indicated in (|43|l . 
In the isotropic limit 7^0 



G+{k) 



V27rmr (1 - i 



if)r(i 



47r3/2 



2 



-ikm 



r 1 



ikm 
m 



r i-i^ 



-i/c In 



k + iQ^ 



-dfc 



l — m 



h^Th 
h:$>Th 



(44) 



In the last line, only the leading behavior due to the 
simple pole at — iO"*" for high fields has been included. 
The integral in equation H44() has been given by Tsvelick 
and Wiegmann^i. It allows to determine the zero- 
temperature scales defined below. There are several sin- 
gularities of the integrand equation in the complex 
fc-plane: 

i) m = /: Poles are distributed in the upper and lower 
half-planes, with an additional dominating cut along 
the negative imaginary axis. 

n) m < I: A cut along the whole imaginary axis goes 
along with sub-leading poles in both the upper and 
lower half planes. 



Singularities in the lower (upper) half plane are relevant 
for h>Th {h< Th). Poles fc„ = i(2ri + 1), n = 0, 1, . . . 
in the upper half plane only give a leading contribution 
in the exactly screened case I = m. They have residuals 
2i( — l)"+^/rt!, such that the magnetization is given by a 
series 



M,{h<Th) ' = ' 



(-ir 



r(n + i) 



^ (2n+ l)n!r(m(n+ i)) 

/ I 1\ \ "'("+5) / , s 2n+l 

Th 



One recognizes the signature of a Fermi liquid in first 
order : Mi[h) cx h and Xi(^ = 0) =const. Upon inserting 
explicit values, one finds agreement with equation l|21|l . 
We shall establish this agreement explicitly for arbitrary 
7 below, equations H55() . H56|l . 

Let us draw our attention to the cut in the lower half 
plane, for values I > m. By linearizing the integrand, we 
find 



M,(/i > Th) 



I i>m ml I 1 mining 
A= i [-m In TO -I- (m- 2) In 2] .(45) 



The contribution ©(In "h/Th) is absorbed by the defi- 
nition 



M,{h > Th) 



Th = 

I l>7 



c~^Th 



(46) 



ml 



In- 



2 In^A 

T,, 



+0 [In-' h/Th 



(47) 



Equation defines the zero-temperature scale for large 
magnetic fields. 

Finally, upon "encircling" the cut in the upper half 
plane, one only replaces the pre-factor / in equation l|47|) 
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by / — m and sets — Inh/Th = hiTh/h: 



m m{l — m) 



ln4r 

^ h 



, In In ■ 



+ O Tnlh] (48) 



Note that this simple replacement does not hold for 7 ^ 
0, as can already be seen from the lowest order, equation 

The free spin value of the magnetization is approached 
logarithmically at high fields. This asymptotical free- 
dom in the weak coupling limit is a genuine feature of the 
Kondo model. An analogous effect occurs for low fields in 
the under-screened case; however, the first correction is 
of opposite sign compared to the high-temperature case, 
cf. equations H47|) . H48|l . Classical Fermi liquid behavior 
appears at low temperatures if the impurity is exactly 
screened. The physical origin of these results has been 
revealed by Nozieresi already before the exact solution 
of the Kondo model was known: At high fields, correc- 
tions to the asymptotical freedom of the impurity spin are 
caused by the weak anti-ferromagnetic coupling with the 
host particles. At low fields, the impurity spin is partially 
screened due to strong anti-ferromagnetic exchange. Two 
kinds of interactions with this impurity-electron system 
may occur: On the one hand, a weak /erromagnetic cou- 
pling of the residual spin S — m/2 with the host, due 
to the Pauli principle (this explains the change of sign 
in the leading order on the rhs of equations H47() . H48|l ). 
On the other hand, a polarization of the bound complex 
by host electrons, analogously to the Fermi liquid excita- 
tions at 5 = m/2. This polarization is given by the poles 
and dominated by the ferromagnetic Kondo interactions, 
reflected by the cuts in the complex plane. 



Finite temperature 

After the calculation of the magnetization, we proceed 
with the specific heat for I > m = 1. By inserting equa- 
tion H40() into equation H33|l . we get 



-—V{-\nh/n), 



3hK,\e[m 

where the Fourier transform of V{x) is given by 

''^''^-\^m[unk) + i], i>i ■ 

From equation H4()(l it follows for I = 1 that V{k) 
G+{k)s{k), so that 



C,{T^TkMIT- 



X{h) , (49) 



where we have inserted equations (|32|l . H39() . (|40|l . The 

constant |e'i(0)| is determined from the condition that 
equation H49(l matches equation H21() for Z = 1 = m, so 
that |e;(0)| = l/{KiG^{0)), and 

C,;(T « Tk, h)/T . (50) 



3.4. Over-screened and exactly screened cases, 
2S < m (ground state) 

In this section, we calculate the ground state for / < m, 
finite temperatures are dealt with in the next section. 
If I < m, there are m equations, the last one being 



A 



where = 2(1-^^/^) • ^^'^ right-hand side of equa- 
tion H28|l . only the m*'' entry is different from zero. Con- 
sequently, 



-e-{k)+d^{k) + ^5{k) 



(51) 
(52) 



Equations (|36|l . H51|l are solved by the Wiener-Hopf 
method, afterwards e+ is inserted in equations (|77|) . ifS^ . 
The relevant matrix entries read 



A. 



[k) 

Al,m{k) 



With A™,„(/c) 
are given by 



2 cosh ^fc sinh sinh ^ ( - - r 

sinh^ sinh 
2coshffc sinh/2^ sinh ^ (7 ~ 



4rXk) 

G+{k) 



sinh sinh ■ , ^ 



G^{k) G-{k), the auxiliary functions 

J_ G+{k)G-{~i) 
2^ (fc + iO+)(/c-t-i) 

(2m^(l-m2))*r(l-i|) 



r(i-i|)r(i 
r(i-if |)c 



-iak 



A gI(% ■ tl^^ of high fields h:^Th, 



and 

the magnetization reads 
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r . r + i|) r (1 + r (1 - i|) r (i - i| .-j^m ^^^^ ^ 

r(i-i|(^-m))r(i-u|)r(i + i^|) ^ + io+ ^'-^2- ^^^^ 



The isotropic limit 7-^0 yields 
G+(fc) 

^ /-^r(i + i|)r(i + imf)r(i-i 



\/2^r(i-i|) 
r(^-i|)r(i-im|: 



ifcm 
~27 



Mi(/i) 



47r3/2 



r(i-i;|)r(i + i/f; 



ikm\ 
~2e"y 



-iA; In -, 



A: + iO+ 



-dA: 



Note that for 7 = 0, a cut along the negative part of 
the imaginary axis occurs. It dominates the poles in the 
lower half of the complex plane, and one is faced with the 
expected Kondo behavior for h > Th, again in accordance 
withi. 

For h <Th, the leading behavior is given by the poles 
with smallest positive imaginary part. 



h \ / h ^ 'Z™-' 



(54) 



Thus over-screening induces non-integer exponents of ft,, 
independent of 7, at T = 0. 

Finally note that for I = to, the two expressions equa- 
tions (|44|l . (|53|l coincide. It is shown that the first non- 
vanishing order linear in h of the magnetization leads to 
the T = 0, /i = susceptibility, calculated in equation 
H21(l . For I — TO, the magnetization reads 



M,{K) 



47^3/2 



^^ r(i + i|)r(i-i|) r(i-iff 
--'r(i-iif) r(i-i-i--m 



-iA; In tjt lak 



-dk 



Being interested in fields h < Th, one takes account of 
the poles at fc„ := i(2n-t- 1) with residuals 2i(— l)"+^/n!. 
These result in the series 



1 



(-1)^ 



2n+l 



7ri/2 (1 - Im) ^ (2n+l)nl 

r(l-t-n)r(^(n+l))a^"+^ 

\(TO(i+n))r((^-TO)(n+i)) 

From the definition of T/j, equation (|30|l . one gets to first 
order in h: 

Tn h 

hm M,{h) = -— — — (55) 



X^{T=Q,h = {)) = 



27r (1 - Im) T, 
1 



K 



211(1- ^to) Tk 



7^ • (56) 



This result is expected from equation H21|l . 



Finite temperature 

As pointed out above, we are not able to account 
for corrections of the linearized functions tj<m, to > 1 
in the framework of the rigorous linearization. In the 
following, we consider the case j < m for small mag- 
netic fields and low temperatures. For 7 = 0, this case 
was treated previously by Affleck^ using bosonization 
and CFT-techniques. He calculated the low-temperature 
Wilson-ratio analytically for arbitrary to. Analogous re- 
sults have been obtained by the TBA-solutionS only for 
TO = 2, 5* = 1/2. Here, we will confirm Affleck's findings 
for TO > 2, extended to the anisotropic case 7 7^ 0. 

First consider the case to = 2, 25 = 1. The relevant 
auxiliary functions are given by: 



In yi{x) 
lnb2(x) 



[s*ln»2 5B2] [x) 

I3h 



^ + [s* fii ^i] (x) + [k * In Sa] (x) 



MT,h) = - — 



[k_ * InSa] (x) 
T_ r°° lnYi(a;) 
27r ./-TO cosh ( X -t- In : 



dx . 



(57) 



In section lSTl we found (equation ipn|l specialized to to 
2): 



limfs* In 052*821 (x) 



Consequently, 



limlnYi(x) = ln2 + e" 



2(7r- 27) 



8 4(7r - 27) 



(58) 



One approximates the kernel in the same scheme as in 
equation 1^01) . 



lim 



2T 



T«T,^ cosh (x + In T/Tk ) Tk 



x<ln^.(59) 



The approximations H58|) . (|59|1 lead to a constant inte- 
grand, which we deal with by introducing a cutoff D, 



MT,h) 
Tk 



T 

An ' 



■In 2 



47r(7r - 27) 



\nT/TK 



Idx , < 
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The cutoff D depends neither on T nor on h, so the oc- 
currence of the logarithm liiT/Tx in the above equation 
leads to 



lim Ci(T,h) 

T,h^TK 

lim Y, (T, h) 



ITIuT/Tk 
4 

1 



InT/Tj 



K 



27r(7r - 27) Tk 



lm=2,2S=l 



3(1 - 27/71) 



(60a) 
(60b) 
(60c) 



confirming, for 7 = 0, numerical findings by 
Sacramento^ and Affleck's prediction. 

For values m > 2, we perform an asymptotic lineariza- 
tion of the NLIE in the region a; — > 00, following^. Cor- 
rections to In Y^.^°°^ are expressed through a correction 
function Dj{x), 

InYAx) = lny/°°' + DAx) , lim DAx) = . 

Linearizing Inj/j to first order in D^, 7 = 1, . . . , to — 1 



lnyj(a;) 
D,{x) 

Dm 



(00)^(00) 



J-1 + 1 
1/2 



i?,(x) 



0111 tt: 



(61) 
(62) 



Equation l|61(l has been derived for a; 0. This 
is approximately accounted for by writing Dj{x) = 
Dj(x)9{x). The linearized equations form an algebraic 
system by Fourier transforming. Each unknown function 
Dj(x) can be expressed in terms of Dm, 

Dj{x) ^[tj*Dm]{x) (63) 
sinh(j + 2)f - sinhjf 



,9(fc).(64) 



As usual, JFfc [<j(a 
tion satisfies equation for j = 1, . . . , m — 2. The 
function g{k) has to be determined from the last equation 
j = TO — 1 and is found to be 



tjik). We see that Dj from cqua- 



2 cos smh TO + 2) — 

TO+2 ^ '2 



(65) 



The first equation, j = 0, is already contained in equation 
(|^ by to = and therefore Dj^q ee 0. Equations 
(I65II imply that the auxiliary functions InY,- behave as 



InYjix ^ 00) - InF, 



(00) 



const. X e 



-4a;/(m+2) 



(66) 



The constant is given through lnS„iQ5m, equations (|62ll . 
(I63|l . and will be determined numerically below. Com- 
bining equations (|63|) . H64|) and (|65|) . one finds for the 



impurity part of the free energy 



MT,h^TK) = -Tln- 



sm — — 

m-t-2 



sm 



m+2 



00 cosh ^ 



^D2sik)e''^''^^/^-dk. (67) 



Since T < Tk, only negative imaginary values for k are 
allowed. 

The leading T-contribution is given by the singularity 
closest to the real axis of the integrand in equation l|67|l . 
For TO > 3, this is a simple pole at kp = —41/ (to -I- 2), so 
that 



MT,h<t:TK) = -Tlii- 



(si 



7r(2S+l) 
m+2 



m+2 



Ta 



T 



47rS • 2ir(S+l) 

sm -^^^ sm — ^ Z. ' 

m-t-2 m+2 



sm sm , „ 

m+2 m+2 



Tk 

47r(S+l 



(m + 2)cos;;^cos;|^sinI^ 

e^ ln[Q5„,»„](x)da; (68) 

The specific heat and susceptibihty are derived: 



C,{T,h = 0) = a\ 



4(m + 6) / T \ 



x.(r,/i = o) = ( ^ 



92 



a(/3/i)2 



(69b) 



h=0 



The derivative in equation Ij69b|) only acts upon the in- 
tegrand in equation H68() . Both C,;/T and Xi show the 
same T-dependence, yielding the low-temperature Wil- 
son ratio, defined in equation H23|l : 



a(r) xhiT) 



4^ 
~3~ 



(to -I- 2y o(f}hy 

4:{m + 6) 



a 



(70) 



/i=0 



for TO > 3. The quantity a defined in H68() depends on 
S, 771, 7 and (3h, which is taken to be zero in (|69a|l . Ij69b|l . 
This means that the a-factors in (|69a|) . (|69b|l are con- 
stant and known, so that from measuring Ci or Xi, the 
scale Tk can be determined. So the parameter Tk has 
the meaning of a low-temperature scale also in the over- 
screened case. 

We did not find an analytical expression of the inte- 
gral in equation (|68|) . but determined the coefficients in 
equations (|69a|l . Ij69bp numerically. Therefore, the NLIE 
are solved as described in the appendix, and the coef- 
ficient of the leading T4/(™+2)-decay of Ci{T < Tk), 
Xi{T <^ Tk) ■ T is extracted from the numerical data. In 
fig. [3 some curves are shown. The Wilson ratio for the 
isotropic case is 



17=0 



(2 + 77i/2)(2 + to) 
18 



■; TO> 2. (71) 
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For finite anisotropy 7^0, the magnetic field is scaled 
by the factor 1/(1 — m-y/'K) in the NLIE, equation (|12|l . 
(|13(l . Furthermore, the non- integer powers in (|69a|l . 
(|69bp are determined by the fusion hierarchy, (|65|l . where 
the anisotropy parameter 7 does not enter explicitly. 
We therefore conjecture that in analogy to the exactly 
screened case (I23|) and two-channel overscrcened case 
(|6Ucf) . the Wilson ratio for finite 7 is 



1 



(2 + m/2)(2 + m)2 



1 — rwy /it 



18 



; m > 2 (72) 



This means that anisotropy is irrelevant in the renormal- 
ization group sense, in accordance with previous treat- 
ments of the underscreened case by TBA2i2ii&. On the 
other hand it has been shown in the CFT-approach24 
that a generic anisotropic spin-exchange is relevant. 
Since the spin-exchange (|2b|l is a polynomial in the spin- 
operators, the situation here is different from that in24. 
We leave further investigations on this issue as a future 
task. In table m our results are compared with equations 

(EH), dzn). 



S=\I2 
S=l 
3=3/1 
S=2 




-40 



-30 



-20 



FIG. 7: The low-temperature Wilon ratios for tfie 
overscreened case with m — 5 and anisotropy 7 = 
0, 0.2, 0.5, 0.87max, where the lowest set of curves is for 7 = 0. 



4. HIGH TEMPERATURE EVALUATION 

In this section, the NLIE arc analyzed in the limit of 
high temperatures T 3> Tk by an asymptotic lineariza- 
tion. The fact that contrary to the TBA approach, we 
deal with a finite set of NLIE, is most useful here. In 
the isotropic single channel case, our results are farther 
reaching than those obtained from the TBA equationai^, 
the treatment of the isotropic multichannel case and of 
the anisotropic models in this regime is novel. 

This section is divided into three parts. In the first 
part, the isotropic limit 7 = is considered. The sec- 
ond part treats the low-temperature limit of the under- 
screened case, which is conceptually very similar to the 
regime of high temperatures. The anisotropic case is 
treated in the last part. 



4.1. Isotropic case 7 = 0, h — 

In the high temperature regime, the first corrections to 
the asymptotic value of the free energy are given by the 
asymptotic corrections of the concerned auxiliary func- 
tions for X —00. Corrections are of algebraic and 
logarithmic-algebraic nature and can be found within an 
asymptotic linearization scheme. This is also true for 
Z > m in the low-temperature limit, referring to the re- 
gion X —^ 00. In this case, the set of I equations de- 
couples into two sets of m and I — m equations and the 
high-temperature-linearization is directly transferable. 

One makes use of the knowledge of the asymptotic be- 
havior of the auxiliary functions derived in appendix^to 
extract the high-temperature behavior of the free energy, 
the specific heat and magnetic susceptibility in zero field. 
Employing once again the approximation equation ljA3|l , 
this time in the integral of the free energy, we find 



lim MT) 

t:»Tk 



lim C,{T) 



lim \i{T) 



, , ml(l + 2)Tr T 

-T\n{l + 1) + \^ ' ^ 3 ^ (73a) 



12 In^^ 



1 ml{l + 2)Tr^ 



(73b) 



l{l + 2) ( ^ m m^lnlnT/Tif 



i2r 



In 



2 T 
Tk 



(73c) 



Values for (j){m) are given in table ^ for m = 1, , 
Let us cite the special case m = 1, 



1) 



0.04707±2-10" 



,5. 



(74) 



Equations H73a|l - H73c|l constitute the first calculation of 
the leading orders of the specific heat and magnetic sus- 
ceptibility for general m in the framework of an exact 
solution. The single-channel case m = 1 agrees with 
known TBA-resultSfi^. Especially, the coefficient of the 
(lnr/TR-)~^-decay of the magnetic susceptibility is de- 
termined, from which we will calculate Wilson's number, 
relating the high-temperature to the low-temperature 
scale. 

From equation H73c|l . we calculate Wilson's number for 
the underscreened and exactly screened cases, I > m. Let 
us first draw our attention to the exactly screened spin- 
1/2 case, I ~ 1 ~ m. WilsonS^ obtained for the zero- 
field spin-1/2 susceptibility by his renormalization group 
approach 



j,im x,(T) 

Tk 
Tk 



1 



0.1032 ±0.0005 



2ttTk Tk 

2-K ■ (0.1032 ± 0.0005) =: 27r^ . (75) 



This ratio relates the low-temperature scale Tk to the 
high-temperature scale Tr-, defined by absorbing the 
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TABLE I: Numerical results for the low temperature Wilson ratios in the over-screened case m > 25'. In the isotropic case, we 
compare our data with AfHeck's CFT approach—, equation I7H . For finite anisotropy, the values given in brackets are those 
from equation 17211 . 



m 


7 = 


[1] 


7 = 0.2 


7 = 0.5 


7 = 0.8 


2 


0.202638 


0.202642 


0.2251582 (0.2251582) 


0.2071899 (0.2071898) 


0.3377373 (0.3377373) 


3 


0.370 


0.369 


0.4102 (0.4104) 


0.490 (0.493) 


0.607 (0.616) 


4 


0.60791 


0.60793 


0.6749 (0.6755) 


0.805 (0.811) 


0.99 (1.01) 


5 


0.928 


0.931 


1.033 (1.034) 


1.23 (1.24) 


1.51 (1.55) 



term 0{x ^) in the asymptotic expansion 



lim Yi(T) 



1 



[ln!BiSi] jx) 
27rT cosh {x + IuT/Tk) ' 



dx 



1 



2T 

uix) 
"1,0 

"1,0 



dj^ [inSiSi] {x^-IhT/Tk) 



Ax) 

"1,2 



(76) 



x = - In T/Tk 



"1,1 



b'f}\n\x\ 



X ■ 



(x) 

1,2 

(x) 

1,1 



x=- \tlT/Tk 



Wilson's number is identified to be 



27r^ = exp — 7^ 



(77) 



Andrei and Lowenstein^ carried out a perturbative ex- 
pansion of the 25* = m = 1 free energy, both for h/T ^ 
and T/h ^ 0. By requiring that in the first case, the re- 
sult should depend on h/Th, in the second case on T/Tk^ 
they deduced the ratio Th/TK- Moreover, they deter- 
mined the ratio Th/Tk from the (conventional) BA. Ar- 
guing that the ratios of the energy scales are universal 
(unlike the scales themselves, which do depend on the 
cutoff-scheme used), they found by combining their two 
results (the analytical expression is due to Hewson^i) 



eC+3/4 
47r3/2 



0.102676. 



(78) 



We generalize Wilson's definition (|77|l to the general spin- 
s' case, in the presence of m channels. A scheme of nu- 
merically solving the integral equations which allows for 
the calculation of the corresponding ratios is given in the 
appendix. For general Z, m, equation H77() reads in the 
notation of section im 



27r^ = exp 



(m) + 



which only depends on m, analogously to the ratio Th/Tj^ 
for T = 0, equations (|45|l . H46|) . This is in contradic- 
tion withS^. There, the Wilson numbers for S arbitrary. 



m = 1 are calculated. The ratio Tk/T^ for T = is 
found by BA techniques for 7 = and agrees with ours, 
equation H3U|) for 7 = 6. Th/Tx is found by conventional 
perturbation theory. The resulting Tk/Tk depends ex- 
ponentially on S{S + 1). We leave this question to be 
clarified. 



By inserting equation H74|l . one gets for 



1: 



0.102678±2-10" 



This result agrees with equations H75|) . (|78|) . In table ITTl 
Wilson numbers for the general I > m case are given. 
Note that for 2^ > m, the susceptibility at low tempera- 
tures can be obtained from that at high temperatures by 
replacing I {I — m). This does not change the value of 
^, which means that in the under-screened cases, only one 
scale (namely Tk) governs the low- and high-temperature 
behavior. 



m 


0(m) 








1 


0.04707 ± 2 


■ 10"^ 


0.102678 ±2-10 




2 


0.0502 ± 1 ■ 


10"^ 


6.4604 ■ 10~^ ±6 


■ 10" 


3 


0.044 ± 1.4 


10"^ 


4.455 • 10"^ ± 8 ■ 


10"" 


4 


0.038 ±1.4 


10"^ 


3.178 ■ 10"^ ± 7 ■ 


lO"'^ 


5 


0.034 ± 1.4 


10"^ 


2.124- 10"^ ±6- 


10"" 



TABLE II: Wilson numbers for / > m. 



The over-screened case to > ^ is obtained from equa- 
tions H73a|) , Ij73 b|) , H73c|l by replacing to I and inserting 
lnYj=i and its derivatives, equations ljA13|l . (|A26() . into 
the definitions of the free energy, specific heat and sus- 
ceptibility. The Wilson numbers are again given in table 
im where now m ^ I, (j){m) — > 4'{l). In analogy to the 
exactly and under-screened cases, Tr- has been identi- 
fied as a low-temperature scale also for the over-screened 
case in equations (|69all . (|69b|l . From the universality of 
the low-temperature Wilson ratio, we expect the Wilson 
number in the over-screened case also to be universal, in 
analogy to the exactly screened case. 
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4.2. Low temperature evaluation in the 
under-screened case 



In this section, we make use of the results of appendix 
13 from which one deduces the low-temperature behav- 
ior in the under-screened case of thefoUowing quantities 
(r«TK): 



X^{T) 



-T\n{l-m+l) 

T m{l — m){l — m + 2)7r 



In' 



3 Tk 
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T m{l ~ m){l - m + 2)ti'^ 



{I - m){l ~ m + 2) 

i2r 



In^ 



m2 In In Tk /T 4m^<j){m) + ir? j \ 



In^^ 



(79) 



The low-temperature behavior for exact screening is 
given in section H^.H where the expected Fermi liquid be- 
havior shows up. Similarly to T = 0, a change of sign 
in the leading corrections to the asymptotic values of the 
susceptibility is observed, cf. equations H73c|l . Its 
physical interpretation has been given in the sequel of 
equation It applies analogously in this case. 



4.3. Anisotropic case, 7 7^ 

According to equation ||HJ), the anisotropy parameter 
is restricted to < 7 < where n = max(Z,TO). Since 
the kernel k{x) decays exponentially in direct space, cor- 



rections to In 25 



(00) 



InF, 



(00) 



J I , ill 1 J are expected to be exponen- 
tially small. Thus it is no longer permitted to replace 
convolutions with s{x) by algebraic multiplications. In- 
stead, let us write equation (jAl|) in direct space, with 
the same notations as equation (|A2|) . however including 
a finite magnetic field from the beginning. First subtract 
the asymptotes (C := Ph/2), 



„ _.smhC(Z + 1) , , 

lim2e —77 — 5i[x) 

x<.o smh CI 

+ [s* Sl-l + K* Si - *5i] (x)} 



hm {— e"^ 



(80) 



Contrary to the isotropic case, the driving term — can- 
not be neglected since all quantities on the rhs of equation 
(I80|l are exponentially small. We did not find a closed so- 
lution to equation H80() . but determine the exponent of 
the leading exponential decay. Therefore first note that 
the equation 



S^^\x) 



-e-" + 



k{x - y)d':^\y)dy 



is directly solvable by Wiencr-Hopf techniques. This so- 
lution relies on the fact that 



sinh^iiit 

2 7 



2cosh^sinh2f 



is factorizable in functions analytical in the upper and 
lower half planes. The leading decay is S^^\x <C 0) ^ 

eir^. We take the solution J^"' as an ansatz for S-, 
where ( is assumed to be small. It is seen that there is 
no further restriction on the leading decay of 6-, such 
that 

(5_(.T < 0) - e^^ . 

However, we did not succeed in determining the coeffi- 
cient. It depends on especially. 



lim[5_+5_](.T) 
/,(T»Ta-) 



Q(r»TK) 

x^(r»TA-) 



(Si + {phfa2)c^'= 
-T Ml + 1) 

+ (5i + {(3hfa2) 



Tk 
T 



l{l 



12T 



1 - "2 T^r 



Tk 
T 



(81) 
(82) 
(83) 



Note that both Ci and Xi show similar decays, contrary 
to the isotropic case. 

If r <C Tk , one argues in close analogy to the isotropic 
case: The decoupling into two independent sets still 
holds. However, the asymptotic values of the auxiliary 
functions for x 00 are related to their a; — > — 00 
counterparts by substituting I ^ (I — m) and scaling 
(3h — > ajSh, a = (1 — ra'y/n)~^, equation (|15|l . This 
scaling affects the susceptibility: 



/,(T«rA-) 

C,(T«Tk) 
X^{T « Tk) 



-T [\tl{1 
(ai 



1)- 



T 

T^. 
(l — m)[l — m 



2)a^ 



12T 



1-02 



Tk 
T 



(84) 



(85) 



(86) 



The constants in equations H81|) - (|83|l differ from those 
in equations l|84|) - (|86|l . For ease of notation, the same 
symbols have been used. However, if we send T — > first 
with h ^ and / > m, the rigorous linearization, section 
Em is done. Then Xi ^ h'^ ln~^ H/Tk, equation l|Ig|l : 
The T-dependent power-like divergence is replaced by a 
logarithmic, /i-dcpendent divergence. 
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5. CONCLUSION 

We presented a novel exact solution to the anisotropic 
multichannel spin-S* Kondo model. The free energy 
contribution of the impurity is given by a finite set of 
(max[TO, 2S']+l)-many NLIE. By analytical and numeri- 
cal studies, we confirm and extend known properties of 
this model. 

The low temperature case is characterized by Fermi 
liquid behavior in the sense that Wilson ratios are de- 
fined. However, only in the exactly screened case I = m, 
the limits T, /i — > commute and Ci{T)/T, Xi approach 
finite values. These values were calculated by the diloga- 
rithm technique and for / > m = 1 by the dressed charge 
formalism, in the framework of a rigorous linearization. 

We analysed the ground state for arbitrary 
anisotropy, spin and channel number and observed 
non-commutativity of the limits ft., T — > for models 
with I ^ m. In the underscreened case, if /i = 0, an 
asymptotic approach to free (/ — m)/2 spin asymptotes 
was recovered for T <^ Tk, formally analogous to the 
T ^ Tfc case, paragraph 14.21 On the other hand, 
performing first the limit T — > while letting h ^ 0, 
a non-integer rest spin occurs for 7 7^ in the under- 
screened case, equation (|43f) . connected to a quantum 
critical pointSiS. 

We analyzed the low-temperature behaviour of the 
over-screened models and found non-integer exponents 
of /i if T = 0, equation (|54|l and of T if ft = 0, equa- 
tions Ij69bp . Ht)9a|) . Especially, we determined numerically 
low-tempcrature Wilson ratios for arbitrary anisotropy, 
confirming and extending results by AfflecliS^. 

At high temperatures, T ^ Tk, the impurity spin ap- 
proaches asymptotically the behavior of a free spin of 
magnitude S = 1/2. Corrections to the asymptotic val- 
ues depend in their amplitude on the channel number m 
and have been calculated analytically. Especially, Wil- 
son numbers relating low- to high-temperature scales are 
determined with the help of a numerical solution. 

We expect that new insight into the multichannel case 
is obtained by generalizing the lattice path integral ap- 
proach proposed inii for the isotropic 5 = 1/2, m = l 
model. To this aim, i?-matrices must be constructed 
which are invariant under the action of gl(2|l), contain- 
ing higher dimensional irreps of su(2)(2S. 

Finally, it should be possible to derive the quantities 
determined numerically in this work also by analytical 
methods. This question is left open for future investiga- 
tions. 



APPENDIX A: ASYMPTOTIC LINEARIZATION 

1. The region .t — > — 00 

We begin with the asymptotic expansion of In 03;, 
ln5B;, / > m in the region x ^00. The case m > I 
is obtained therefrom afterwards. Consider the equation 



for in the ft = 0-case: 

Inb; = s * InF/-! + K * InS/ - K- * In^i (Al) 

We shall show that ln*Bi(x), lnS/(x) approach their 
asymptotes x — > —00 as x~^ and calculate the corre- 
sponding coefficient 63 ^ 

The first I — 1 integral equations determine the j- 
depcndence of lnyj (a;), j = 1, . . . ,1 — 1. To see this, 
define 

lnY^{x) = 2ln{j + 1) + 2Sj{x) 
ln^i{x) = ln{l + 1) + 5i{x) 
lim 6j{x) = 

limln2/,(x) = lnj{j + 2) + 2^-^^^S,{x)+0{6'^) 

limlnbzfx) = lnl+^-^Si(x)+0(Sf) . (A2) 

ln5B; (Inb;) is related to ln*B; (Inb;) by complex conju- 
gation. The crucial approximation is 

[s*S,]{x)^S,{x)/2 + 0{5'J{x)). (A3) 

Since s is an exponentially decaying kernel, this approx- 
imation is justified for algebraically decaying Sj. Such 
an algebraic behavior is indeed expected from the inte- 
gration kernel k{x), which itself decays algebraically for 
7 = 0. In the anisotropic case 7^0, this is not true, 
since Sj{x) also decays exponentially. The Sj satisfy the 
recurrence relations 

2 S, {x) = 6,+,{x) + (x) (A4a) 

So{x) = (A4b) 

'-^Si{x) = Si^iix) + 3{l + l)d{x) , (A4c) 

where in the last line the function d{x) is to be deter- 
mined and its prefactor has been chosen by convenience. 
These equations determine Sj up to a constant factor. 
Note that from equations (j A4ap . HA4bp 

Sji^) - j(j + 2)d(x) . 

Summarizing, 

InYjix) = 21n(j + l) + 2j(j + 2)rf(a;) 

ln'8/(2;) = + I) + 1(1 + l)d{x) . (A5) 

The functions *B;, 25; are related by complex conjuga- 
tion, the Yi are real-valued. 



In «B; =: Bi,i + iSs,; , In S; =: Bij - 182,1 ■ (A6) 

Define the sum and the difference of the integration ker- 
nels. 
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Asymptotically, 



k'-"^ I X- + 1- 



1 Itt 
2^ ~ 2^ 



oo 



(A7) 



The convolutions with the K-kernels are written in the 
following way: 



K * In 05; — K_ * In 05; = k'l'-' * Sj.i + i k'l'' * i?i,2 



is) 



{x) = k'-'^-' ( X ± i— ) , u ~ s,d . 



The first non- vanishing term in an asymptotic expansion 
of Bi i(x) around \x\ — s- cx) is 



Bi^i{x) ^ - m + 1) + In 



/ + 1 



/ — rn + 1 



ei-x) 



The (^-function has to be understood asymptotically for 
large x. This regime is equivalent to small fc- values in 
Fourier-space, Tk[0{—x)] = —i/k + TrS{k). In this re- 
gion around the origin in Fourier space, 7rA:J-/c[K^'*']/2 ~ 
Thus it follows that 



(x) - In ■ 



/ + 1 



I -m+1 V 4a;3 ix^ 



It is useful to define correction terms to the asymptotic 
behavior for a; ^ 0: 



"1,3 



. "2,2 



lim \n^i{x) = ln(/ + l) + ^+i^ . (A9) 
One then performs the asymptotic expansion 
[x) 

/oo poo \ 

Bi,2{x) dx + b^^U K,(x)dx)(A10) 

/oo poo \ 

Bi^2ix)dx + b2^2 / Ksix)dxj , 

where 

K^'\x)dx = 1 

? + l 



i /I 

TT /I 



Bi,2{x) dx 



-(to ± 10-^) + In 



Z - 771 + 1 



The last integral is done numerically with the indicated 
precision. In the following, we set 

£B„(.)d. .f(-,„ + ln^^) (All) 



Insert equation ljA9p in equation ljA2p and keep only the 
linear order in 5i, 



lnbKx)^ln/ + iii P^ + i^ 
I \ x-^ x^ 



(A12) 



Combining equations (|A8|) . (|A10p . ljA12p one expands 
equation HA1|I around x — > — 00: 

/ + 1 6^'' 

In/ + — ^ ^ = In/ + (/ - 1)(/ + l)d{x) 



I 



4 111 , -^(- r B2 i{x)dx + b^^\) 

4x3 /-m+1 x^\2j_^ ' 



/_+! 6^ 

/ x2 

1 

+ 2 



TT, / + 1 

— In 

4 I -m + 



B2j{x)dx + 62,2 
Using equation (|Alip . we find 



d(x) 



^1,3 



71" TO (;) 

12 x3 ' - 
/(/ 4- 2)m TT^ 
3 T 



lim lny,(x) = 21n(j + l)+j(j + 2) 

X — > — 00 

lim lnS;(x) 



6x3 ' 



(A13) 



HI + 1) - ^-rzT ^ ^ — 



4x^ 



12x3 



Note that the x-dependence of the corrections is deter- 
mined through the asymptotic behavior of the kernel in 
the convolutions K*ln05;, K*ln05;. The amplitudes fol- 
low from the y-hierarchy. 

We proceed with the asymptotic evaluation of 
df3h lnyj (x), dfjh ln*B;(x) in the regime x ^ 0. The InY, 
are symmetric with respect to Ph: The system of NLIE 
remains the same upon replacing ph — > —^h and sub- 
stituting 05/ by 05; . So djSh In Yj vanishes identically for 

h ~ 0. Define B;™' := d/sh In 05;. Then the only equation 
which remains is 



B 



(m) 



(X) = ( 



= I 1 - c"'''®'^^^ 



2 + * 



-K_ * B 



(x) 



(A14) 



'"^ and i?!™^ are related by change of sign and complex 
conjugation. 



B 



(m) 



iB. 



(m) 
1,2 



Below it is shown that the imaginary part of van- 
ishes for /i = 0, so Bi"^^ = B;™'* is real valued. We shall 

determine asymptotic corrections to up to the or- 

der O (x~2), so that corrections to ln05/(±cx)) (of order 
O (x~3^)^ can safely be neglected. In this approximation. 



lim fi_c-'"'»'(^) 



/ + 



t(^ + ^4^) 
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lim fi-c- 

v—>-oc \ 



I 



l(i + ii^)(Ai5) 



I — m + 

and one finds the asymptotic equation for i?^ ™^ , 



^(1 +«'•..«;•■■'). (AM) 



The imaginary contributions vanish as expected. An 
equation analogous to (|A16|) does not exist in the TBA- 
approach, where one deals with the infinitely many Yj 

and their derivatives. But Y^*-™^ = OVj, as mentioned 
above. 

The constant asymptotic behavior of b\™'^ is written 
in the compact form 



b["'\x) = A+ - A_sgna; 



(A17) 



A+ + A_ 



A+- A_ 



I — m 



Equation (|A17|I and similar equations in the following 



have to be understood asymptotically, for 
similar manner, 



oo. In a 



is') 

k} ' * sgnx 



sgnx 

X 

sgncc In I 



*(2) 



1 

X 



In \x\ 



In \x\ 



(Al8a) 
(A18b) 



'^{x) is the digamma function; ^^{2) = -C + 1, C = 
0.577 ... is Eulcr's constant. The asymptotic evaluation 
of the convolutions is done by using distributions. To 
determine accurately the a;~^-coefiicient, the precise be- 
havior of the auxiliary functions around the origin must 
be known. This is done numerically as shown in the ap- 
pendix. 

We make the following ansatz for , which consists 
in extrapolating the asymptotic behavior over the whole 
axis with the aid of distributions: 



(x) = 



A+ - A_sgna; -I- -(A^-P 

X 

A-i- — A_sgna; 1 



AL^^sgnx) -I- i^(A« - A«sgnx) + \{A^:' 



A_^sgna;) 



— A_ 
2x 



A^^^ - aL^^ sgnx 



-(aL')-a(/' 



.(i) \ Inbl 
A_ sgnx 



2x2 



aL^^^'(2) 



A^^' - AL^'sgnx + I 



1 

2^ 



(A19) 



(A20) 



where e is determined numerically (cf. appendix, equa- 
tion ljB3|) ). The convolutions in the second equation are 
done with the help of equations (|A18ap . l)A18bp . The 
terms sgna; and sgnx/x give the next higher order contri- 
bution when convoluted with the kernel. The asymptotic 
form of equation ljA16p is found by inserting equations 
comparing coefficients, one finds 



aV'^ + aL' 



a«-aL' 



(0 



A 



+ 

(2) 
+ 

(2) 



A^ 



(l - m)A_ 



m{l-m) ^(1) 



-lA 



(1) 



■ I 



^ ^ - 8 
= ?(AL^'*(2)+e/2) 

= (Z-m)(AL^^*(2) + e/2) 



(A21) 
(A22) 



In the high-temperature regime, the x ^ — IhT/Tk "C 



behavior is of importance, 

(„) I ml m'^l lii\x\ 

lim B, (x) = — — 

' ^ ' 2 Ax 8 x^ 

I fm^- *(2) e\ 

Compare equation ljA23p with the low temperature re- 
sults equations (|47|l . H48|) . Both are formally identical up 
to the order In Ixj/x^, x = — IuT/Tk in equation ljA23p . 
X = -Inh/fh in equations lI37|), lO- 

We shall focus on the x +00 asymptotes in the next 
paragraph. 

Our analysis is continued by expanding dp/^lnYj =: 

Yj''\ dj^ln^i =: b\^\ a|^ln»i =: B^i'\ These func- 
tions are given by the system 



y/'^) = (1 - e- '"^-^ ) s * [Yj^l + Yj^l 



Y 



(x) 



B 



(x) 



In 



B, 



I 



(A24a) 
(A24b) 

+ (1-e- 



+ (1 _ e- 1" ) * ^(x) _ ^_ ^ ;g(x)^ (^24c) 
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(x) 



■.B 



ix) 
i.i 



(x) 

1,2 



2. The region a; ^ oo 



The vanishing of Bj™'* helps to expand b[^'' : 



(x) 
i.i 



B 



(m) 
i.i 



I 



2{l + l) 



Y, 



(x) 



B, 



l + l 
(x) ^ Trm 
4x2 



Re ( kL"^ * Sj^' 



i.2 



(m) 



Z + 1 8x2 



(x) 
(x) 



/ + 1 



Im ( kL^) * 



(A25a) 



(A25b) 



From equation (|A25b(l . = O (x"^). Together with 
equation (|A7|I , one concludes that the last term in brack- 
ets in equation l|A25a|l is O {x~^) . This means that for 
our purposes, the convolutions in equation (|A25a|l can be 
entirely neglected. Thus in the asymptotic limit, equa- 
tions l|A24a|) - (|A24cll are simplified considerably: 



(X) _ j(j+2) 



Yr 



(X) _ 



B 



(x) 



2(J 


1 r 
7 



+ 1)2 V + J + 1 



B 



(m) 



I 



2{l + l) 



-Y 



(x) 



This system bears similarity with equations (jA4a| 
From the solution of those equations, we conclude 



lim Y^^\x) 



-*(2) 



jU + 2) 
6 

9 

TO 



1 



TO to2 In |x| 

X 



2e 



2x2 
(A26) 



l{l + 2 



lim bW(x)=&^„^ + ^ 



X ''2 ^2 



^(X) 



12 



1 



TO 

X 



-*(2) 



m" In |x| 
2x2 



+ 2e 



In the appendix, a procedure for numerically determin- 
ing the coefficient A^' — A'^^'sgnx of the x^2_£iecay is 
described. It deviates slightly from the analytical esti- 
mate. We introduce an extra symbol (j), 



(m) := 



A 



(2) 



A 



(2) 



= i.f^vi.(2) + £ ) (A27) 



Numerically, is found to be independent of /. However, 
our data do not suffice to exclude a dependence on to. 
The results are given in tabled in the main body of this 
work. 

If / < TO, one replaces in the whole preceding analysis of 
this appendix / ^ to and m —f I. The thermodynamical 
quantities are given by Yi and its derivatives. 



We determine the behavior of the auxiliary functions 
for large values .x — * oo in the under-screened case 
I > m. In this limit, the set of NLIE decouples in 
two separate sets with exponential accuracy. Consider 
first the field-free case h = 0. The last / — m func- 
tions, namely Y^+i , . . . , 25; , *B; , satisfy a set of equations 
formally identical to equations (|A4ap - (|A4c|l . with j, I re- 
placed by J — TO, Z — TO. The whole analysis of the preced- 
ing paragraph n applies to this case; of special interest 
are now the corrections to the x 
of the auxiliary functions. The results are 



In ^ ^ 0-behavior 



<8z(x) = ln(Z-TO + l) -i 



.m{l — to) tt 



4x2 

m{l — to) (l — m + 2)7r2 



lny,>,,„(.x) = ln(/ + l-j) + 



12x3 

m{l — m){l — m + 2)n^ 



6 x^ 



Bt\x) 



f/^^(x) 



bI''\x) 



I — m m {I — m) (l — m) In |x| 

2 ^ 4^ 8 x2~ 

I - m f m2vI/(2) e 
~8 ^ 2 

(j — ?7l)(j — TO + 2) / 771 TO2ln|x| 

6 V 



2x2 



(/ — m)(/ — TO + 2) / 771 77i2ln|x| 

12 V ^ ~ 2^ 



APPENDIX B: DETAILS OF THE NUMERICAL 
TREATMENT 

From the analytical analysis of the NLIE, numerically 
ill-conditioned terms were found, namely those C(l), 
0{x-^) in the limit | x| — > oo^ Since in Fourier space, 
they would appear as simple poles and discontinuities in 
the origin, they are subtracted from the auxiliary func- 
tions before performing the FFT and treated separately. 
After having applied the inverse FFT, the analytically 
convoluted terms are re-added: 

ft * InF = K * (InF - f) + K * / , (Bl) 

N A 

lim Iny(x) - /(x) , (B2) 

\x\ — *oo 

and K is a kernel, In Y an auxiliary function. The function 
/ is chosen to be asymptotically equal to Iny, equation 
(|B2|) and contains the numerically ill-conditioned contri- 
butions. The difference InK — / is therefore numerically 
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transformable. The term denoted by N is treated nu- 
merically, that labeled by A analytically. The Fourier- 
transforms of / and k are known. The convolution fc * / 
is solved separately, such that no further manipulations 
in Fourier space are necessary. We solve numerically 
for (Iny — /): By including terms 0(1), 0{x^^) and 
0(ln l^l/a;^) in f{x), the function (InF — /) decays as 
a;"^. Since the kernels also decay as x~'^, the coefficient 
of the x"^ decay of Iny is obtained from the integral of 
(InF— /), a quantity which is accessible numerically with 
high accuracy. Consider the case In y = B^"^^ , given by 



(IA14|I . and define the regular function Breg ■= B^' — /, 
such that Bj.eg{\x\ — > oo) ^ x~'^. Then for |x| — > oo 

/oo 
Breg{x) dx 
-oo 

/oo 
k{x) dx , (B3) 
-oo 

and the quantity e in equation (|A20|) is identified as e := 

J^^Bregix)dx. 
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